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Abstract. We apply non-extensive statistics, namely Tsallis statistics, on a case of 
supercritical string cosmology (SSC) as the one studied in pQ and obtain interesting cosmological 
modifications. At the beginning, we seek for the non-extensive corrections to the dilaton energy 
'■/""i ' density and the off-shell (off-critical) terms and by using the set of dynamical equations of [T], 

we derive the modified evolution equation for the radiation energy density in a r.d.e. (radiation 
dominated era). This in fact is characterised by fractal( "exotic") scaling and this seems to be 
a generic result of our analysis. The modified Boltzmann equation was also considered , giving 
us the effects on dark matter candidates relic abundances of non-extensitivity besides the SSC 

■ effects. An effort to give a physical interpretation to the essential results of models such as the 
Q_) ' Tsallis statistical model has been done in collaboration with N. Mavromatos and S.Sarkar by 

■ using the D-particles foam model which is presented at the last part of this work. 



C — , 1. Introduction:Tsallis statistics framework 
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Tsallis statistics formalism is based on considering non-extensive (non-additive) entropies, which 
can reproduce the standard Boltzmann-Gibbs entropy at a certain limit of the free parameter of 
the theory,known as the non-extensive parameter, often denoted as q.The important thing about 
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these families of entropies is that, when applied to ordinary statistical mechanics, they give rise 
to probabilities that follow power laws instead of the exponential laws of the standard case( for 
details on this see [2])- In fact this last property makes Tsallis approach a very convenient tool in 
the attempt to describe systems with long range interactions, e.g gravitation, or systems having 
long memory effects. In most cases that Tsallis formalism is adopted ,e.g.[3], the non-extensive 
parameter is taken to be constant and very close to the value for which ordinary statistical 
mechanics is obtained (q = 1), meaning that all calculations can be performed in the leading 
order to (q — 1), with no considerable loss of information. These two assumptions are also 
made throughout our work whereas in general it would be interesting to consider, at some point, 
a time-dependent parameter q, with a function dependancy that could be treated as a fitting 
parameter of our model. 

It is also important to mention that although Tsallis entropy gives probabilities which have 
physically, a different behaviour than in the standard case, it is constructed in such a way 
that it preserves all the "good" properties of the Boltzmann-Gibbs entropy. Namely, it is 
always positive, it is concave which is crucial for thermodynamical stability and it also preserves 
the Legendre transform structure of thermodynamics. This in fact is what differentiates Tsallis 
approach with other attempts made to walk away from extensitivity. 



2. Tsallis statistics effects on SSC: a case study 

We consider the SSC model adopted in [T] and seek for the modifications due to non-extensitivity. 
The set of dynamical equations in the Einstein frame in this case, will have the form: 
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3H 2 =p m +p^ + ~—G^ (1) 

2H = -p m — f>(t> — Pm — Pcf) + (2) 

<p + 3H(f) + + -(p m -3p m ) = ---^ - —G,p (3) 

where all densities now carry the non-extensive corrections. We also make clear that in the 
non-critical case, the central charge deficit that comes into the potential for the dilaton field is 
not constant but evolves with time through the Curci-Paffuti equation (for details see [I] and 
[1]). The equilibrium energy densities for matter (non-relativistic limit) and radiation are given 

in [3]. By using these expressions, as well as a correction of the form ef rdt where T accounts for 
the off-shell, dilaton terms and can be overall seen as a source term, we can explicity write the 
contributions of the ordinary matter to the energy density p m appearing in (JT|) , d2J) , (EJ) . 
For instance the non-relativistic energy density will be given by: 



t 



pUt) = ^mm^ie-^l^ + + { ^nl™ (4) 

where it is easy to recognize the standard expression and the two corrections. 

In equations (H|) , ([2]) and ([3]) we have neglected any non-extensive corrections to the off-critical 
terms and Ga, which is a valid approximation under the assumption that those terms are of 
order less than {q — 1). Therefore all that we are left with and always to leading order in {q — 1) is 
to determine the non-extensive corrections to the dilaton energy density and the exotic matter 
which is also taken into account in the density p m by the authors of pQ . By assuming a radiation 
dominated era and by introducing a convenient effective number of degrees of freedom defined 
by the relation: 

9eff, q = 9eff, q + ^ T~*Ap (5) 

where Ap incorporates all the dilaton and off-shell terms, which are assumed not to be 
thermalized and by using for g e ff,q the expression given in [3]: 



9eff,q = ^2 9i,bosons(7^) +g ^2 9j,fermions(^) +(#-l)[9.58 ^ 5i,6osons(^r) 4 +8.98 gj 

(6) 

we can obtain the correction to the dilaton field energy density: 



AM = 4> 2 + + fjro iriV ~ 1)[9-197J2 9bosons + 8.623 9fermions] (7) 
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For the exotic matter we considered that any q-dependance will come into its equation of state 
parameter that anyway is treated as a fitting parameter in the numerical analysis performed to 
©,©,©( see also [1]). 

With all these in hands we can write down the modified continuity equations stemming from 
combination of (P),© and ©: 



^ + 3H(p m + p m ) - <P(p m - 3p m ) +6HC + ^ = 6(H + 0)% 
at E dt E a 

where by C we have denoted the q-correction to the dilaton energy densityi.e: 
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Now from ([8]) we can easily obtain the evolution equation for radiation which will have the 
form: 
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dt E 3HH E Ht 2 E 3Ht E 3iJ2 

where the parameter e appearing above is defined through: 

p r 

One can try to solve (jlOh perturbatively in e since it contains (q — 1) which can always be 
made arbitrarily small by chosing appropriately the non-extensive parameter q. We obtain: 
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where we have set: 



p r {t E ) PS [pr,in a 4 ] e = pr,in a 
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Prom (I10p and the definition for 5 it is clear that under the assumption of non-extensitivity 
and for the case that we are studying, radiation scales in a fractal way and this result arises 
naturally in our analysis. 




Figure 1. : Plot of the radiation energy density xa 4 with respect to the Einstein time 



In the above plots we have used the numerical estimation e ~ 0.2885(g — 1) calculated for a 
r.d.e. 



2.1. N on- extensive effects on Boltzmann equation 

In all the following analysis, we restrict our attention to supersymmetric dark matter species, 
such as neutralinos, which are thaught to be the lightest supersymmetric particles (LSP). In fact 
such an analysis can be related to a quite interesting and rich phenomenology (see also [5] ) and as 
such it can be a very promising tool. In order to describe the number density evolution of those 
species, one needs to write down the modified Boltzmann equation due to non-extensitivity and 
the off-shell/dilaton terms (we recall that our case study is a SSC case). Those last corrections 
for the standard case, are given in [5]. 

The corrected Boltzmann equation for a particular species i/j of mass and in terms of the 
parameters Y^ = — (number per entropy density) and x — nas the form: 

dY I VH~ ^ 

— = 0.264 (a A \v\)gj f h q m m pl [(Y^ q ) 2 - (Y^ eq ) 2 } Y^ q (14) 

where h q denotes the entropy degrees of freedom and its form in the non-extensive case is given 
i" ®- 

By looking at (|14p we see that this has exactly the same form as in [5] with the only difference 
that the number densities as well as the effective number of degrees of freedom and the entropy 
degrees of freedom all carry the non-extensive corrections. We note here that the last term in 
(114ft represents the contributions from the extra terms (off-shell and dilaton) which are elegantly 
incorporated in the T term. 

Following the standard approach in solving the Boltzman equation, let us assume that the 
species under consideration decouples from the plasma at a specific temperature Xf (freeze-out 
point). It is clear that at temperatures greater than Xf, the interactions of the species maintain 
it in thermal equilibrium, therefore Y^ tq = Y^ q e q at all those times. Substituting this in (|14p . 
yields: 
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Y^ q ( X ) = 0.U5g s (h q )- 1 x-2e-x[l + ^-(- + 3x- 1 +X~ 2 )]e (15) 

and g s counts the particle's spin degrees of freedom. 

In order to determine the freeze-out point, we use the general freeze-out criterion: T(xf) ~ 
H{xf) and we also evaluate 1151 at Xf- Therefore we obtain the "new" freeze-out point that will 
be given by: 

Xf,q = Xf~ \^Xf l +X7 2 ][-T + 3 Xj 1 +Xj 2 ~ -r-^— (9.58^2 g ijbosons + 8.98^2 g jJermions )} 

11 4 ' 9eff,q 

_ (16) 

where all effective number of degress of freedom are counted at the freeze-out point. We also 
note that by Xf m <fTSl we mean the "standard" freeze-out point satisfying the relation ([5]): 
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Xf 1 = /n[0.038 {a A \v\) Xf g e / fXf ~ 2 m m pl g s ] + -Zn(^L) + J (17) 

From (|1T|) we see that the non-extensive correction to the freeze-out point of the species if 
can be either positive or negative depending on how early or late the species decoupled from 
the plasma. Very roughly we could say that if the species decoupled at early times (large Xf) 
where also the relativistic contributions were quite large, then with the non-extensitivity taken 
into account the species should have decoupled even earlier, whereas for small freeze-out points 
(late decoupling) the species in the non-extensive case should have decoupled even later. 



With all these in hand and by following exactly the same approximations as in [5], we find 
that today's relic abundances for the species ip will be given by the expression: 
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where the first term stands for the standard result: 
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with J = J {a a \ v\) dx and we have defined as in [5], 

XO 



-JLMZl dx 

Hx) = Xe xo X (20) 



The last term is the correction braught by the non-extensitivity and has the form: 
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f - 1 i ( 1 1 u „ „ \ W H) (x/, g - x/) . x/, g - x/ n ? rg- 1 

/c ™ ti °- - 1+( ^^"^7) )(5e// ' 9 " 5e//) j + ^x7T (1+ i w" x) 
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whereas the second and the third term in (|18|) represent the corrections due to the off- 
shell/dilaton terms estimated by the authors in [5]. 

To do a rough estimation: the modification braught by the SSC term (third term in (|18p ) 

is of the order of — ~ 10~ 13 if matter contribution is considered to be negligeable and we also 

x/ 

take the off-shell terms to be very small (see [5]) whereas in the last term of (]18p the biggest 
contribution comes with a factor (q — 1) which in general is greater than 10 -13 . Therefore 
the effect of non-extensitivity on today's relic abundances can be significant for appropriate 
choices of q and as such can be used to put constraints on the model by comparison with the 
corresponding phenomenology. 



3. D-particles foam model 

In this part we briefly present a physical model that as we show, can give rise to cosmological 
results similar (not necessarily identical) to those of Tsallis statistics. The basic idea for this 
model, is that the spacetime is viewed as a "fuzzy" entity, revealing this way its dynamical aspect 
imposed by a general relativistic point of view .The origin of this "fuzziness" within the context 
of quantum gravity, is found in the interactions of point-like stringy defects of spacetime (D- 
particles) with closed strings. Then as a result of the momentum transfer during the interaction, 
there is a D-particles recoil and this can happen in a stochastic way, meaning that any fraction 
of the momentum of the string could be transferred to the D-particles. Due to this process, 
the spacetime behaves as a randomly fluctuating environment and this is expressed through 
microscopic fluctuations of the metric. Therefore if one considers the case of flat spacetime, a 
generic metric would have the form: 

/ -1 km k 2 r 2 k 3 r 3 \ 
km 1 
9flu ~ k 2 r 2 1 

V k 3 r 3 1 / 



(22) 



where in this model we have chosen the off-diagonal metric components g° l to depend on the 
momentum of the particle as well as on a set of statistical parameters (ri,r2,rs) that describe 
the amount of momentum that is transferred in each collision to each spatial direction. It is 
important to note that these parameters are assumed to follow a gaussian (normal) distribution 
with a standard deviation <r, that is not constant but is distributed with a x 2 distribution , i.e. 
its average value will be: 



<<Ji>= J <Jif(<Ji)d<Ji = a i0 (23) 
o 

with the probability density given by: 

/(<ri) = f(i)W 2<r - (24) 

In order to derive the statistical description for particle systems in our model, we need to 
obtain an expression for the particles energy in a quantum level approach. Let us take the 
massive Klein-Gordon equation in a general gravitational background: 

[g^d^dy - m 2 ]$ = (25) 

If we expand this form and we also assume plane wave solutions, we obtain the following 
expression for the particle energy: 



2g 0i ki ± ^A(g 0l ki) 2 - Ag 00 (g u {h) 2 + m 2 ) 



^2 = ( 26 ) 



where we have thrown away all cross terms, that is terms of the form kikj (see also [6]). 

Subsequently, we find what (|26p gives for the case of the metric (|22|) . The only approximation 
we make is to drop again all cross terms. In order also to simplify our expression we make 
the substitution: E 2 = k\ + k\ + k 2 + m 2 and this in fact would correspond to the particle 
energy in the case of a Minkowski spacetime. But in our case ( (|22[) ) this is just a NOTATION 
and nothing more since the actual particle energy is given by the expression (|26p . Now with 
the calculated energy in hand and by using the standard definition of statistical mechanics: 
< n >= — rm 1 — vm, we can obtain the distribution functions for fermions and bosons in our 
model: 



< n > 



exp(f3(E - n)) + £ 
exptfjE-n)) T(f + 2) 2 
(expi^E-^ + O 2 T(§) n 2 

[(P 2 + ^)(a 2 m kt + a 2 2 kl + a 2 3 kl) + ^m 2 (a 2 01 k 2 + a 2 2 k 2 + a 2 3 k 2 )] (27) 

where £ = +1 applies to fermions and £ = — 1 applies to bosons. By recalling that aot has 
dimensions of inverse temperature (see [22]) it is convenient to set: aoi = where aoi is now 
dimensionless. Under this substitution, ([27]) becomes: 

1 

< n > 



exp(/3(£-/x))+£ 
exp(^- M )) T(§ + 2) 2 
(exp((3(E-ri) + a) 2 r(f) n 2 

[(/3 4 + ^)(a i 2 kf + a 02 2 kl + ah 2 4) + ^m 2 {a m 2 k 2 + a 02 2 k 2 + a 03 2 k 2 )} (28) 



This is our expression for the distribution functions from which one can reproduce all number 
and energy denisties for bosons and fermions in both the relativistic and non-relativistic limit. 

4. Conclusions 

We saw that when Tsallis statistics was applied to one case of SSC we obtained quite interesting 
cosmoligical effects. At first we obtained an "exotic" scaling for the radiation energy density 
since our analysis referred to a r.d.e, which seemed to be a rather generic effect. In fact if the 
same analysis had been performed for a m.d.e, fractal scaling for "dust" matter should occur 
as well. The second interesting effect has been the modifications to the relic abundances of 
CDM candidates (neutralinos) which seem to dominate on the effects due to off-shell/dilaton 
terms presented in [5]. These modifications should in fact be constrained by the relevant 
phenomenology but this is a part of another work. In the last part of our analysis, we saw 
that by using this modified metric (I22p we obtained a distribution function expression ( (|28p 
that is very similar to the result obtained in [3], where the modified distribution functions were 
derived by introducing a non-additive parameter q to the standard Boltzman-Gibbs entropy 
(Tsallis entropy). The smallness of the deviation from the standard case (1st term in (|28p ) was 
guaranteed in our case by the smalness of the statistical parameters coi whereas in the Tsallis 
statistics case this came from requiring small departure of q from the value 1. Therefore, we have 
proven that we managed to get a physical result very similar to the result of Tsallis statistics by 
simply starting from a different background spacetime physically stemming from stringy defects 
(D-particles) recoil. The same analogy holds for all the other quantities defined in[3] since we are 
performing our analysis with equivalent distribution functions. However those results will appear 
in a future work together with the effects on dark matter relic abundances since in this part of 
our analysis our main goal was restricted to demonstrating the analogy between the physical 
effects of Tsallis statistics and the D-particles foam model. 
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